Frustrated order on extrinsic geometries 
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We study, analytically and theoretically, defects in a nematically-ordered surface that couple to 
the extrinsic geometry of a surface. Though the intrinsic geometry tends to confine topological 
defects to regions of large Gaussian curvature, extrinsic couplings tend to orient the nematic in the 
local direction of maximum or minimum bending. This additional frustration is unavoidable and 
most important on surfaces of negative Gaussian curvature, where it leads to a complex ground 
state thermodynamics. We show, in contradistinction to the well-known effects of intrinsic geom- 
etry, that extrinsic curvature expels disclinations from the region of maximum curvature above a 
critical coupling threshold. On catenoids lacking an "inside-outside" symmetry, defects are expelled 
altogether. 
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The intrinsic curvature of a surface frustrates the or- 
der of materials living upon it, making defects necessary 
even in the ground state This frustration derives 
from the incompatibility of straight and parallel direc- 
tions on surfaces with Gaussian curvature. A nematic 
texture on a sphere, for example, must have a net topo- 
logical charge of +2 0-3], as do, for example, lines of 
latitude on the globe. On more complex surfaces with 
non-uniform intrinsic curvature, the ground-state orga- 
nization of defects is well-described theoretically by com- 
petition between defect-defect interactions, which favor 
separation of like-charged defects, and curvature-defect 
interactions, which favor the localization of defects in re- 
gions of relatively high Gaussian curvature [H, 0] ■ This 
picture is far from complete, however, as real surfaces 
are endowed with an extrinsic geometry that depends in 
detail on how they sit in space [7H9( . Stripes on a curved 
surface, for example, tend to orient along "flat" directions 
to minimize their ben ding energy, an effect observed in 
block copolymer films [10| . 

In this Letter, we demonstrate that the extrinsic geom- 
etry of a surface generates a strong and previously un- 
studied source of frustration which controls the ground 
state of anisotropically-ordered surfaces. While the in- 
timate relationship between topological defects and in- 
trinsic curvature is well-established, we show that the 
coupling between the spatial organization of defects and 
extrinsic geometry is as important, and in some cases 
more important, than intrinsic geometry for determining 
the ground state of nematic order on anisotropic sur- 
faces. Away from umbilics, extrinsic geometry gives rise 
to symmetry-breaking fields in regions of high curvature 
that favor uniform order most strongly in precisely those 
regions where the intrinsic geometry favors the locally 
"isotropic" textures generated by defects. The frustra- 
tion between textures favorable to intrinsic and extrinsic 
energies becomes particularly pronounced on surfaces of 



negative Gaussian curvature. As prototypical example 
this frustration, we study the ground state of nematic 
order on the surface of catenoid (see Fig. rjj . While the 
intrinsic curvature favors the confinement of four —1/2 
disclinations to the neck of the catenoid, the extrinsic 
curvature favors local orientation along the directions of 
zero or maximum curvature along the catenoid, which 
tends to repel defects from the neck. As the relative 
strength of extrinsic to intrinsic effects is increased, this 
competition results in either a discontinuous or contin- 
uous expelling of the defects away from the neck, and 
in some cases off of the surface completely. Anisotropic 
ordering on surfaces of complex curvature is common 
to many materials systems-from bicontinuous particle- 




FIG. 1. (color online) (a) A catenoid with four, partially ex- 
pelled — 1/2 defects as determined by numerical minimization. 
Defects are confined to the neck in the ground state of intrin- 
sic energy, shown in side- and top- view in (b) and (c). In 
the ground state of extrinsic energy, (d) and (e), defects are 
expelled from the necked results in an overall +1 charge on 
the planar regions. 
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stabilized emulsions, or "bijels" 14j . to biological mem- 
branes enriched by curvature-inducing proteins [l5j|-and 
we expect the frustration between extrinsic and intrinsic 
effects to play a generic role in shaping the structure and 
thermodynamics of these systems. 

We begin by briefly overviewing the effect of intrin- 
sic geometry on curved nematics. We define an orthog- 
onal coordinate system on the surface u = (u, v) and 
a corresponding orthonormal frame {e„,e„,N}, where 
N is the unit surface normal. A two-dimensional ne- 
matic phase, therefore, has a director lying in the tan- 
gent plane such that n = cos#(u)e„ + sin0(u)e„. For 
convenience, we choose a local isothermal, or conformal, 
coordinate system, (it, v), defined by a metric of the form 
ds 2 = n(u)(dv 2 + du 2 ). In the presence of a configura- 
tion of disclinations, labeled by n and having topological 
charges s n and positions (v n ,u n ), one finds the explicit 
expression 6(u) = — ^ rl (s„/27r)Imln(u + iv — u n ~ iv n ) 
and arrives at an energy written in terms of disclination 
position alone [f|, 0] , 
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(1) 



where C is a constant modulus equal to the Frank con- 
stants in the one-constant approximation. The first term 
describes long-range interactions between disclinations, 
where Vi n t (u) is given by inverting the Laplace-Beltrami 
operator, V 2 K:„t(u) = S 2 (u). The second term de- 
scribes the interaction of the disclination with the Gaus- 
sian curvature, G(u) , as encoded by a "geometric poten- 
tial," V G (u), satisfying V 2 V G (u) = G(u) fl. This high- 
lights the fundamental consequence of intrinsic geometry: 
disclinations are attracted to regions of oppositely-signed 
Gaussian curvature via long-range interactions. 

Thus, defects, through the intrinsic energy, "screen" 
Gaussian curvature, which would otherwise force lines 
of parallel orientation to converge or diverge. This is 
neatly demonstrated on the catenoid, a surface which 
concentrates negative Gaussian curvature within a nar- 
row "neck" adjoining two asymptotically planar sur- 
faces. Conveniently, the radial and azimuthal direc- 
tions provide a natural isothermal coordinate system, 
respectively, with V 2 = [I/f2(u)](<9 2 + d 2 ) and 
G = -V 2 (lnO)/2. We immediately obtain V G (u) = 
-Tnf2(u)/2. For a catenoid with neck radius, r, 
we obtain the geometric potential explicitly, Vc(u) = 
— In [cosh 2 (w/r)] /2. The intrinsic energy thus favors 
negatively charged disclinations along the line of max- 
imal curvature, v — 0, and positively-charged defects are 
expelled. The geometric potential is strongly confining 
for negative disclinations as Vg ~ — \v\/r for \v\ 3> r, 
indicating the divergent cost of expelling these disclina- 
tions to infinity. Intuitively, the strong confinement of 



these defects can be understood from the uniform 9 con- 
figurations shown in Fig. Q] (b-c), which due to the ab- 
sence of the requisite balance of negatively-charged de- 
fects in the neck effectively carry +1 disclination tex- 
tures on each of the adjoining asymptotic planes. Be- 
cause J dA G = — 4-7T, a net defect charge of —2, or four 
— 1/2 disclinations, are required to screen Gaussian cur- 
vature of the neck. 

Extrinsic geometry further complicates this already 
complex picture. The extrinsic energy contributions to 
the free energy quantify the cost of out-of-plane gradients 
of n and are, therefore, sensitive to the surface's embed- 
ding, X(u). The surface curvature can be expressed as 
a tensor, hij = N • c^c^X, where the indices can take on 
the values s or u. The eigenvectors of the tensor give 
the directions of maximal and minimal curvature while 
the corresponding eigenvalues give the principal curva- 
tures. In the most general form, the coupling between 
and the director can be expressed in terms of the 
three rotationally invariant quanitites: C nn = tfnPhij, 
C nt =n*(Nx n) j hij and C tt = (N x n)*(N x n) ] h l3 [n|. 
To second order in curvature, all terms allowed for achiral 
materials include [7[ 

E ex = J dA{K nn C 2 n + 2K nt C 2 t + K tt C 2 t + 

+ K' nn C nn + K' tt C tt }. (2) 

For a microscopically symmetric interface such as a 
tilted-bilayer membrane, for which N — > — N is a sym- 
metry, terms linear in the curvature must vanish (i.e. 
K' nn = K' tt = 0). For asymmetric interfaces, however, 
linear terms account for the physical distinction between 
"inside" and "outside" Unlike the intrinsic energy, 

which is invariant under global rotations of the director 
around N, the extrinsic coupling generates an unavoid- 
able and geometrically-induced symmetry breaking field 
on anisotropic surfaces. In terms of the angle field, 



dA{j cos[4(6> + p)] + i cos[2((9 + /?)] } (3) 



The strengths of the 4- fold and 2-fold symmetry-breaking 
fields are determined by the local anisotropy of bending 
on the surface, 7 = (K u + K nn — 2K nt )(ni — K2) 2 /8 and 
and i = (K' nn - K' tt )( Kl - k 2 )/2 + (K nn - K tt )( Kl - 
K2)H. Here, the principal curvatures of the surface are 
Ki and K2, H = («i + K-x)j1 is the mean curvature, and 
(3 measures the angle between e u and the direction of 
largest principal curvature. For the case of the catenoid, 
note that both H = and f3 = 0. 

From this point of view, extrinsic energy leads to a 
generic preference to lock the director to the direction 
of local maximum or minimum curvature, a preference 
which is strongest in regions of high-curvature anisotropy 
where n\ — K2 is large. For example, the curvature dis- 
tribution on the catenoid, k\ = —K2 = r~ 1 sech 2 (?j/r), 
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leads to the extrinsic preference for uniform orientation 
most concentrated in the neck region. Importantly, the 
uniformly-ordered ground states of the extrinsic energy 
on negatively-curved surfaces like the catenoid [Fig. Q] 
(d-e)] correspond to the maximal energy configurations 
of the intrinsic energy. 

To explore the defect configurations that best negoti- 
ate the compromise between intrinsic and extrinsic ener- 
gies on the catenoid, we employ two methods. First, we 
calculate the perturbative contributions from E ex in the 
limit cases of small 7 and 7'. To first order, this amounts 
to evaluating Eq. ([3]) at the intrinsic energy saddle-point 
with an additional global rotation of 6 determined by in- 
trinsic energy minimization for 4 disclinations distributed 
at intervals of 7r/2 around the neck and staggered at 
heights v — ±h away from the neck. Second, we numer- 
ically study the ground state of a coarse-grained model 
of nematics on the catenoid. Our approach is based on a 
generalization of the Lebwohl-Lasher model [l3| for the 
case of a two-dimensional mesh of non-uniform geometry 

Starting with a catenoid discretized in a large num- 
ber of small triangular patches, a nematic director is as- 
signed to the tangent plane of each patch. The spin-spin 
interaction energy is given by H int = J2i,j AAiAA.j [l - 
(rvn,-) 2 ] V(ry), where V(rij) = exp[-r?-/ (2a 2 )} weights 
spin coupling between patches i and j by their spatial 
distance ry, AAi is the area of the i th patch and we 
choose er = 0.36r. The area weighting has the important 
advantages that the coarse-grained model appropriately 
reduces to the continuum form of the intrinsic energy 
in the limit that A A; 0, and it minimizes the influ- 
ence of nonuniform mesh geometry. To this we add a 
finite-difference approximation of Eq. §3§ to account for 
extrinsic energies. A Monte- Carlo simulated annealing 
scheme with the Metropolis-Hastings sampling method 
is used to determine the ground state 

For symmetric interfaces, 7' = 0, the relative strength 
of the extrinsic and intrinsic energy scales is character- 
ized by AK/C, where AK = (K nn + K tt - 2K nt )/8. In 
Fig. HJa), we plot the mean vertical distance measured 
from the neck at v = along the symmetry axis of the 
surface of the four -1/2 disclinations in the ground state. 
Both the simulations and perturbation theory show that 
the defects are confined to z — by the geometric poten- 
tial for small extrinsic couplings, |Aif| < AK C 0.8C. 
At this critical coupling, the four disclinations are found 
to jump discontinuously to configurations outside of the 
high-curvature neck. We denote this as a transition be- 
tween the confined and partially expelled configurations, 
as the defects in the latter state, though expelled from 
the neck remain bound by the geometric potential and 
maintain a finite distance from the neck for any finite 
value of AK [Fig. HI (a)]. 

For an asymmetric interface with 7 = 0, we find dra- 
matically different behavior. Here, the relative strength 
of extrinsic to intrinsic energies depends on scale, char- 




AK'r/C 



FIG. 2. (color online) Mean vertical distance of the defects 
from the neck of the catenoid for a (a) symmetric and (b) 
asymmetric interface. The solid curves are perturbation the- 
ory predictions, and the points are numerical simulation re- 
sults. 



acterized by the ratio rAK'/C, where AK' = (K' nn — 
K' tt )/2. As shown in Fig. HJb), instead of a discontinuous 
transition, the defects are pushed from the neck region 
for any finite value of AK' , with mean distance growing 
linearly as r\AK'\/C for weak extrinsic couplings. Most 
significant, the four defects are expelled to v — > ±00 as 
|Aif'| — > AK' C ~ C ' jr and are completely absent from the 
ground state of nematic order for stronger anisotropics. 
In this case, we see a transition from a partially expelled 
to fully expelled state for the disclinations. 

These two results for symmetric and asymmetric in- 
terfaces demonstrate the profound influence of extrinsic 
geometry on the ground state of nematic order. They 
also highlight, in different ways, the subtle relationship 
between disclinations and the anisotropic extrinsic cur- 
vature. In the simplest view, the extrinsic anisotropy is 
localized to a finite region of size r in the catenoid and, 
therefore, competes with the near-field isotropic texture 
that surrounds a defect core. A disclination in the neck 
costs an extrinsic energy of roughly r 2 AK(An) 2 m AK 
relative to the uniformly-oriented state. Relocating the 
disclination away from the neck adds an intrinsic energy 
penalty of order sC. Comparing, we expect defects to 
be expelled when AK 3> C. However, J dA(An) 2 = 
— J dAG = 47T, so there is only a finite extrinsic energy 
to be gained; disclinations can only be expelled partially. 

This simple picture is not adequate on asymmetric in- 
terfaces. Indeed, for nematics on asymmetric catenary 
interfaces, the extrinsic director coupling is far more ho- 
mogeneous in its distribution. This is straightforward 
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FIG. 3. (color online) Conformal map of defects on a catenoid 
a distance h from the center of the neck (v = 0), shaded 
by extrinsic energy density with streamlines tangent to the 
nematic director. Integration contours C v < h and C v < h 
(dashed lines) have clockwise orientation. 



where p is the radial distance from the axis and z(p) is the 
vertical height from the neck. Since n(p) ~ d 2 p z ~ p~ a ~ 2 
and the area element grows as dA » pd(f)dp, it is straight- 
forward to estimate the total extrinsic energy gain for 
expelling defects on asymmetric surfaces is, roughly, 
/ dAn ~ J fl dpp~^ a+1 \ where R is the system size. For 
a > this cost is finite and, when compared to the intrin- 
sic energy cost to expel disclinations of ln(R/r), suggests 
that defects are always confined on the surface. On the 
other hand, the cost diverges when R/r — > oo for a < 
suggesting that defects are unstable for any finite extrin- 
sic coupling. Interestingly, the catenoid manifests the 
marginal case of a — 0, in which the partially-expelled 
and fully-expelled states are separated by a second-order 
critical point. These simple arguments suggest that the 
ground-state thermodynamics of defects on curved inter- 
faces is critically sensitive not only to the microscopic 
physics underlying the extrinsic and intrinsic couplings 
but also to subtle changes in surface geometry. 



to see in isothermal coordinates, where the area ele- 
ment grows exponentially with vertical height from the 
neck as dA = dudv cosh 2 (v/r). This precisely balances 
the |Ak| = r _1 sech. 2 (i;/r) decay, indicating that extrin- 
sic coupling is strong everywhere. We can understand 
the far-field interactions between defects and extrinsic 
coupling by visualizing a partially-expelled nematic tex- 
ture, with defects at v = ±/i, on the conformal map 
(the v — u plane) (see Fig. [3|). The intrinsic energy 
is AK'r^ 1 J dudv cos[20(u)]. We consider contributions 
from constant-?; contours by examining the rotation of 9 
around closed contours, such as C v> h and C v< h shown 
in Fig. [3] . In the map, the director rotates by +2w rel- 
ative to the upper and lower boundary to maintain the 
uniform order in the asymptotic planes. The value of 
§ c di ■ X70 — 2ir — irn(v), where n(v) is the number of 
disclinations in the region v' > v. Since the contribu- 
tions from the vertical portions of the contour cancel due 
to the u — y u -\- 2nr symmetry, the director rotates by 
+2tt along contours for \v\ > h, leading to no net extrin- 
sic energy gain. For contours \v\ < h, no net rotation of 
9 occurs, leading to a uniformly-oriented texture in the 
region spanning the disclinations and a coherent extrin- 
sic energy gain of roughly —AK'h. Comparing this to 
the intrinsic cost of pulling the defects from the neck, 
Ch/r, wc find that the extrinsic effects are deconfining 
when AK' > C/r, and disclinations are expelled entirely 
from the surface. Notice that the size scale of the surface 
becomes critically important for determining the ground 
state on asymmetric catenary surfaces. 

We conclude by briefly noting that the competition 
between extrinsic and intrinsic effects is generic. For ex- 
ample, we may consider the class radially-symmetric sur- 
faces with the topology of the catenoid whose geometry 
far from the neck is described by z(p 3> r) = r(r/p) a , 



We thank R. Selinger for computational resources and 
fruitful discussions. This work was supported in part 
by an allocation of computing time from the Ohio Su- 
percomputer Center. This work was supported as part 
of Polymer-Based Materials for Harvesting Solar Energy, 
an Energy Frontier Research Center funded by the U.S. 
Department of Energy, Office of Science, Office of Basic 
Energy Sciences under Award Number de-sc0001087 



[2 
[3 
[4 

[5, 
[6 

[7] 
[8] 

[9 

[10 

[11 

[12 

[13 
[14 



badel@polysci.umass.edu 



grason@pse.umass.edu 
x csantang@physics.umass.edu 
[1] V. Vitelli and A. M. Turner, Phys. Rev. Lett 93, 215301 
(2004). 

T. C. Lubensky and J. Prost, J. Phys. II 2, 371 (1992). 
D. R. Nelson, Nano Letters 2, 1125 (2002). 
H. Shin, M. J. Bowick and X. Xing, Phys. Rev. Lett. 
101, 037802 (2008). 

D. R. Nelson and L. Peliti, J. Phys. 48, 1085 (1987). 
V. Vitelli and D. R. Nelson, Phys. Rev. E. 70, 051105 
(2004). 

W. Helfrich and J. Prost, Phys. Rev. A 38, 3065 (1988). 
P. Biscari and E. M. Terentjev, Phys. Rev. E 73, 051706 
(2006). 

H. Jiang, G. Huber, R. A. Pelcovits and T. R. Powers, 
Phys. Rev. E 76, 031908 (2007). 

C. D. Santangelo, V. Vitelli, R. D. Kamien and D. R. 
Nelson, Phys. Rev. Lett. 99, 017801 (2007). 
These invariants can be constructed from the other two 
and invariants of the curvature tensor, H and G. 

E. P. Lewandowski, J. A. Bernate, P. C. Searson and K. 
J. Stebe, Langmuir 24, 9302 (2008). 
P. A. Lebwohl and G. Lasher, Phys. Rev. A 6, 426 (1972). 
E. Kim, K. Stratford, R. Adhikare and M. E. Cates, 
Langmuir 24, 6549 (2008). 

[15] J.B. Fournier, Phys. Rev. Lett. 76, 44364439 (1996). 



